Probing Noise in Flux Qubits via Macroscopic Resonant Tunneling 
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Macroscopic resonant tunneling between the two lowest lying states of a bistable RF-SQUID is 
used to characterize noise in a flux qubit. Measurements of the incoherent decay rate as a function of 
flux bias revealed a Gaussian shaped profile that is not peaked at the resonance point, but is shifted to 
a bias at which the initial well is higher than the target well. The r.m.s. amplitude of the noise, which 
is proportional to the decoherence rate I/T2*, was observed to be weakly dependent on temperature 
below 70 mK. Analysis of these results indicates that the dominant source of low frequency (1//) 
flux noise in this device is a quantum mechanical environment in thermal equilibrium. 

PACS numbers: 85.25.Dq, 03.67.Lx 



1 

-a 
c 

o 
o 



> 
OO 

m 
00 
O 

o 



X 



The viability of any scalable quantum computing archi- 
tecture is highly dependent upon its performance in the 
presence of noise. In the case of superconducting qubits 
it has been shown that low frequency (1//) flux noise 
is of particular concern [1]. Furthermore, evidence sug- 
gests that such devices gcncrically couple to an ensemble 
of effective 2-level systems (TLS) that may be materi- 
als defects [2] . A number of theories exist that attempt 
to correlate these two observations [3], however it is not 
certain whether the TLS observed in spectroscopy exper- 
iments are the dominant source of low frequency noise in 
these devices [4]. The development of additional exper- 
imental probes of low frequency noise will prove critical 
in the quest to build reliable superconductor-based quan- 
tum computing hardware. In this article we demonstrate 
a new experimental procedure for quantifying low fre- 
quency flux noise in RF-SQUID qubits. The procedure 
developed herein complements other approaches to using 
qubits as spectrometers for studying noise [5]. 

Macroscopic resonant tunneling (MRT) [6, 7] is an im- 
portant probe of quantum effects in Josephson junction- 
based devices. In an MRT experiment flux tunnels be- 
tween two wells of a double well potential when energy 
levels are aligned. The tunneling rate and width of the 
tunneling region are strongly influenced by the environ- 
ment. The effect of flux noise on the two lowest energy 
levels can be described using an effective Hamiltonian 
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where e is the bias energy between the wells, A is the 
tunneling amplitude, Q z is an operator that acts on the 
low frequency modes of the environment and cr x r z -\ are 
Pauli matrices. In general, a transverse coupling to the 
environment should also exist but it is believed to be 
subdominant to the longitudinal coupling in flux qubits 
[8] . It is demonstrated that the experimental results pre- 
sented herein are consistent with this expectation. 



A theoretical analysis of MRT in the presence of low 
frequency (non-Markovian) flux noise was reported in 
Ref. [9]. The transition rate Toi from state |0) to state 
|1) (both eigenfunctions of a z with eigenvalues =pl, re- 
spectively) was found to be 
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where S(uj) = ± j dt e iuJt (Q z (t)Q z {0)) is the unsym- 
metrized noise spectral density [10]. Tio is obtained by 
substituting e p — > — e p . Equation (2) is valid provided 
the noise spectrum is peaked at low frequency and the 
integrals defining e p and W are finite [9] . The latter con- 
straint may require low and high frequency cutoffs. 

The quantities e p and W represent the energy shift 
and width of a Gaussian shaped tunneling rate, respec- 
tively. W is the r.m.s. amplitude of the noise and e p 
is a measure of the asymmetry of S(oj). For a classical 
noise source S(u))?ztS(— oS), hence e p = 0. For a quan- 
tum source S(u>) need not be symmetric and therefore 
e p 7^ 0. In thermal equilibrium at temperature T e g the 
fluctuation-dissipation theorem dictates that S(u>) can 
be expressed as a sum of symmetric and antisymmet- 
ric components S(u>) = S s (lu) + S a (to), where S s (u>) = 
S a (w)coth(w/2T off ) (h/k B = 1). If S a {uj) is sharply 
peaked near ui = then it can be shown that 



du>S a (u>)(2T cS /uj) = 2T eS e p 
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The width W is an important parameter in adiabatic 
quantum computation [11] as it defines the precision to 
which a target Hamiltonian can be specified in the logical 
basis defined by the eigenstates of a z . From the perspec- 
tive of gate model quantum computation W is closely 
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related to the decoherence time T 2 * in the energy basis de- 
fined by the eigenstates of Eq. (1). For a low frequency 
environment the decay due to dephasing has the form 
e -t-/2T 2 ^ w h crc l/Tjj = cos (77) W and r) = arctan(A/e). 

To demonstrate the MKT method for characterizing 
low frequency noise we employed an RF-SQUID qubit 
[12, 18]. Previous experimental observations of MRT in 
RF-SQUIDs have been limited to tunneling into higher 
energy levels with and without the help of microwave ac- 
tivation [7, 13]. In this paper we examine MRT between 
the two lowest energy levels of an RF-SQUID. Measure- 
ments of directional tunneling rates r i and r 10 are used 
to extract quantitative information regarding S(u)). 

A schematic of a compound Josephson junction (CJJ) 
RF-SQUID is shown in the inset of Fig. 1(a). It con- 
sists of a main loop and CJJ loop subjected to external 
flux biases <&% and §>5P , respectively. The CJJ loop is 
interrupted by two nominally identical Josephson junc- 
tions connected in parallel with total capacitance C q and 
critical current I q . The CJJ and main loop possess in- 
ductances L c ^ and L q , respectively. If L° n <C L q then 
the RF-SQUID Hamiltonian can be written as 
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(b) Qubit Potential 
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where $ q represents the total flux threading the main 
loop, IF is the conjugate momentum, Ej = ^oI q /2TT 
and $0 = h/2e. This device can be operated as a qubit 
for e [0.5, l]$o and <&% f=a 0. Denoting the ground 
and first excited state of H r f at $| = by |<?) and 
|e), respectively, the qubit states can be expressed as 
10} = (Iff) + |e» Isft and |1) = (| 5 ) - |e» The bias 

energy of Eq. (1) is given by e = 2 \l q \ <S>|, where the 
persistent current \I q \ = |(0| $V i9 1 0> | = |(1| 3> q /L q |1)|. 
The tunneling amplitude of Eq. (1) is given by A = 
(e|W r f|e) — (g\Hrf\g). Both \l q \ and A are controlled 
by S^jp. Maximum A ~ ui p , where ui p is the plasma fre- 
quency of the RF-SQUID, is obtained at = $ /2. 
For $^jp w $0 one expects A — > and the system be- 
comes localized in |0) or |1) . In this latter regime \I q \ 
generates a measurable amount of flux that can be de- 
tected via an inductively coupled DC-SQUID read-out 
(not shown) as described in Ref. [15, 16]. 

The CJJ RF-SQUID from which the data presented 
herein were obtained was fabricated on an oxidized Si 
wafer using a Nb trilayer process with wiring layers iso- 
lated by sputtered SiC>2. The device parameters were 
L q = 661±6pH, C q = 146±3fFand/ c 9 = 1.95±0.05 fxA. 
L q was measured using a breakout structure and C q was 
inferred from the value of L q and measurements of the 
MRT peak spacing [7]. From the MRT spacing we also 
determined oj p ~ 10 GHz. The DC-SQUID was observed 



FIG. 1: (Color online) (a) MRT flux control bias sequence as 
a function of time. Inset depicts a CJJ RF-SQUID subjected 
to the external biases and "3>J. Sense of the macroscopic 
persistent current states |0) (counterclockwise) and |1) (clock- 
wise) are noted, (b) Evolution of the qubit potential U(<& 9 ). 



to have a maximum switching current 7^ = 1.9 ±0.1 /zA 
and the readout-qubit mutual inductance was M ro . q = 
16.7±0.2 pH. The device was mounted in an Al box in a 
dilution refrigerator and all on-chip cross couplings were 
calibrated in-situ as described in Ref. [16]. 

Our experimental procedure is a variant of the MRT 
technique first developed by Rouse, Han, and Lukens [7]. 
We exploit <3?^p to modulate A using a bias line with 
bandwidth ~ 5 MHz. A depiction of the control sequence 
and the evolution of the qubit potential [/($ 9 ) arc shown 
in Fig. 1. The state of the qubit is initialized by setting 
<§> c v w $o/2 (i) and then tilting [/($«) via to an ini- 
tial value of $1 (ii). Slowly raising $^ to $0 in the 
presence of the tilt traps the system in its groundstatc 
(hi). With tunneling suppressed, & q is reset to a target 
value <&j! (iv). Thereafter is lowered (v) to Q C P for a 
prescribed amount of time t during which the system can 
tunnel from the initial state to the lowest lying state in 
the opposite well (vi). Raising $£p (vii) to $0 then local- 
izes the qubit state in |0) or |1) (viii) which can then be 
distinguished by a single shot readout. The probability of 
tunneling from the initial to final state is then measured 
as a function of t and & q . For a given the probability 
of the system being found in |0) can be calculated from 
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FIG. 2: (Color online) (a) MRT from the initial state |0) ver- 
sus t at T t h = 28mK. Results are shown for <3>^ = — 0.488 m<E>o 
(solid squares), — 0.014 m$o (hollow circles) and 0.554 m"I>o 
(solid circles). Slopes of the linear fits yield — roi($'). (b) 
Toi (hollow) and r w (solid) versus $« at T th = 28mK. The 
peaks nearest = have been fit to Eq. (2). 



balancing T i and r 10 : dP /dt = -T Q1 P a (t) + T w P 1 (t), 
where Po(t) + Pi(t) = 1. In the limit t — » the system 
starts in a definite state and this expression reduces to 
dPo/dt = -T i (Po(0) = 1) or r 10 (Pi(0) = 1). 

The results shown herein were generated using a value 
of <&y° for which we measured = 0.56 ± 0.02 /zA. 
This particular target C J J bias was chosen as 1 /T varies 
by nearly four orders of magnitude (10 fj,s— » 100 ms) as a 
function of in the vicinity of = 0, which then takes 
full advantage of the dynamic range of our apparatus. 
Using the calibrated device parameters and Eq. (4) it 
was determined that the above mentioned value of \l q \ 



could be achieved for $ J 1 



0.606±0.001$ . The tunnel 



splitting was then estimated to be Ao = 0.10_g'57 mK. 
The asymmetric error bars on Ao are a consequence of 
its exponential sensitivity to errors in L q , C q , and 1%. 

We have measured P Q ($%.,t) for $j E [-3,3] m$ and 
for the qubit initialized with -Po(0) = 1 and -Pi(O) = 1. 
Example decay data are shown in Fig. 2(a) for Pq(0) = 
1 at three nominal target values of $1 and the device 



FIG. 3: (Color online) (a) Temperature dependence of e p 
(solid) and W (hollow). W has been fit to a line (dotted) 
and e p to a power law C/T& (dashed), (b) T cB = W 2 /2e p 
versus Tth- Dashed line indicates T e g = T t h- 



thermalized at Tth = 28 mK. A summary of initial decay 
rates as a function of <3?^ for both initializations is shown 
in Fig. 2(b). We have independently calibrated = 
by measuring the population statistics in the limit t—^oo 
and observing that the center of the resultant thermal 
distribution is independent of initial state. Fitting of 
the distribution to Po(£^oo)=i _ tanh (e/2T t h)] yields 
T t h- Note that Toi and Tio arc mirror images across 
<&j! = but both are asymmetric about this point. The 
data suggest that Toi and Tio consist of broad peaks that 
are displaced away from $ ^ = in the direction opposite 
that in which the system was initialized. These peaks can 
be attributed to splitting by noise of a single MRT peak 
corresponding to tunneling between the two lowest lying 
states of the bistable RF-SQUID. This splitting indicates 
that S a (u)) 7^ 0. Beyond |<&j|>2m$ , the decay rates 
increase due to MRT between the initialized state and 
the first excited state in the opposing well. 

The data in Fig. 2(b) have been fit to Eq. (2) using W, 



e p and A as free parameters. The fact that the tunneling 
rate can be fit to a Gaussian lincshapc indicates that the 
noise is dominated by low frequency components. We re- 
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lowest energy states of a bistable RF-SQUID have been 
used to characterize flux noise. Analysis indicates that 
the noise source is a quantum mechanical environment 
in thermal equilibrium whose spectral density is sharply 
peaked at low frequency. 

We thank J. Hilton, G. Rose, C.J.S. Truncik, 
A. Tcaciuc and F. Cioata for useful discussions. Samples 
were fabricated by the Microelectronics Laboratory of 
the Jet Propulsion Laboratory, operated by the Califor- 
nia Institute of Technology under a contract with NASA. 



10 20 30 40 SO 60 70 

T th (rnK) 

FIG. 4: Temperature dependence of A. Dashed line is a fit. 

peated the measurements for five different temperatures 
and the resultant fit values of e p and W are summarized 
in Fig. 3(a). Here it can be seen that W is only weakly 
dependent upon T . Fitting these results to a line indi- 
cates that W(T — * 0) w 80 mK. This conclusion corrob- 
orates the observations presented in Ref. [17] from phase 
qubits in which the Rabi decay time is observed to have 
a weak T dependence as T — > 0. In contrast, e p varies 
strongly with T, behaving as - i/T t ° h 80±0 05 w 1/T. Fi- 
nally, an explicit demonstration of the self-consistency of 
our analysis is shown in Fig. 3(b) where T e g obtained 
from Eq. (3) is plotted versus T t h- Agreement between 
T e fi and Tth demonstrates the equilibrium nature of the 
low frequency flux noise in our system. 

Figure 4 shows the variation of A with temperature. 
Note that the results are comparable to Ao as estimated 
from qubit parameters, but that the uncertainty in Ao 
makes it difficult to draw any quantitative conclusions 
from such a comparison. Nonetheless, the data suggest 
that A oc Tth- It is known that high frequency flux noise 
can lead to renormalization of A [10] and that a poly- 
nomial dependence of A on T results if the environment 
possesses an ohmic density of states at high frequency 
[18]. A similar combination of sharply peaked low fre- 
quency noise, as required for Eq. (2), and ohmic high 
frequency noise has been reported for other qubits [1, 19]. 

In addition to studying MRT between the two low- 
est lying states, we have measured MRT from the lowest 
state in the initial well to higher levels in the final well. It 
was observed that the tunneling rate versus could be 
fit to a lincshape whose width increased monotonically 
with the number of levels below the target level, as an- 
ticipated in Refs. [6] and [9]. The width of the peak for 
tunneling to the 50-th level in the final well was roughly 
twice that for the lowest order MRT peak. This weak de- 
pendence on the number of states below the target level 
indicates the relative weakness of S(u>) at high frequency. 

Conclusions: Measurements of MRT between the two 
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